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We  present  a  full  band  structure  calculation  of  temperature-  and  wavelength-dependent  two-photon 
absorption  (TPA)  coefficients  and  free  carrier  absorption  (FCA)  cross  sections  in  InAs,  InSb,  and 
Hg|_tCdvTe  alloys.  The  wavelength  dependence  of  the  TPA  coefficients  agrees  well  with  a  widely 
used  analytical  expression.  However,  the  magnitudes  of  the  TPA  coefficients  obtained  here  are 
smaller  by  a  factor  of  1.2-2. 5  than  the  analytical  values.  In  addition,  the  TPA  coefficient  is  found  to 
depend  sensitively  on  the  photoexcited  carrier  density  in  small  gap  materials.  The  FCA  is  found  to 
arise  predominantly  from  hole  absorption.  The  FCA  cross  section  is  found  to  be  independent  of  the 
carrier  density,  but  is  strongly  dependent  on  the  temperature.  The  calculated  TPA,  FCA  coefficients, 
and  lifetimes  are  fitted  to  closed-form  expressions  and  are  used  in  solving  the  rate  equations  to 
obtain  the  transmitted  pump  and  probe  intensities  as  functions  of  incident  intensity  and  sample 
thickness.  The  calculated  pump  transmission  and  time-dependent  probe  transmission  in  InAs  agree 
very  well  with  the  measured  values.  ©  2007  American  Institute  of  Physics. 

[DOl:  10.1063/1.2718874] 


INTRODUCTION 

The  nonlinear  optical  absorption  and  refraction  proper¬ 
ties  affect  the  propagation  of  high  intensity  light  through 
transparent  materials.1'”  Several  physical  mechanisms,  in¬ 
cluding  two-photon  absorption,  free-carrier  absorption,  and 
thermal  and  free  carrier  refractions,  are  responsible  for  non¬ 
linear  properties.3  6  Both  theory  and  experiments  are  con¬ 
stantly  being  improved  to  accurately  evaluate  nonlinear  ab¬ 
sorption  and  refraction  coefficients. 7~-2  When  the  energy  of 
high  intensity  incident  photons  is  more  than  half  of  the  band 
gap  energy  of  the  semiconductor,  the  electrons  in  the  valence 
band  (VB)  absorb  two  photons  to  reach  the  unoccupied  states 
in  the  conduction  band  (CB).  Then  the  photoexcited  elec¬ 
trons  and  holes  further  absorb  a  photon  with  or  without  the 
assistance  of  a  phonon,  causing  further  reduction  in  intensity 
as  the  light  travels  through  the  material.  Since  the  two- 
photon  absorption  (TPA)  plays  a  vital  role  in  initiating  the 
nonlinearity  in  the  material,  a  number  of  studies  have  been 
undertaken  to  accurately  evaluate  the  TPA  coefficient.  How¬ 
ever,  the  extraction  of  TPA  values  from  these  experiments 
often  depended  on  the  ability  to  control  other  factors  that 

simultaneously  affected  the  nonlinear  absorption.  Conse- 

7—10  15  18—21 

quently,  the  extracted  values  varied  substantially.  ’  ’ 

On  the  theoretical  side,  Wherrett’s  expression  “  with  nonpa¬ 
rabolic  correction11  was  used  widely.  Later  it  was  modified, 
and  the  interaction  matrix  element  was  parametrized  based 
on  the  ability  of  the  kp  models  to  successfully  predict  the 
trends  in  various  semiconductors. 14-1 7  This  expression 
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served  an  excellent  role  in  providing  qualitative  trends  and 
giving  approximate  values  for  the  TPA  coefficients  used  in 
understanding  light  propagation.  However,  an  accurate 
evaluation,  with  full  band  structures  and  matrix  elements  in- 
eluded,  of  TPA  coefficient  made  recently  has  revealed 
that  the  calculated  values  differ  substantially,  by  as  much  as 
a  factor  of  5,  from  those  given  by  the  simple  analytical  for¬ 
mula. 

Light  propagation  is  further  affected  by  free  carrier  ab¬ 
sorption  (FCA),  in  which  the  photogenerated  photoexcited 
carriers  absorb  photons.  The  cross  section  for  the  free  carrier 
absorption,  defined  as  the  ratio  of  free  carrier  absorption  co¬ 
efficient  to  the  free  carrier  density,  is  generally  assumed  to  be 
a  constant  in  light  propagation  studies.3  6  However,  in  our 
detailed  calculations  we  found  that  the  cross  section  is 
strongly  dependent  on  wavelength  and  temperature.  We  have 
carried  out  these  calculations  with  full  band  structures  for 
InSb,  InAs,  and  Hg , _  tCd  vTc  alloys.  The  band  gap  energies  of 
InAs  and  InSb  at  300  K  are  0.35  and  0.175  eV,  respectively. 

In  this  paper,  we  first  provide  a  brief  description  of  the 
band  structures,  the  quasi-Fermi  energies,  and  the  nonlinear 
absorption  coefficients  (TPA  and  FCA).  The  wavelength  and 
carrier  density  dependences  of  the  TPA  and  FCA  coefficients 
in  InAs,  InSb,  and  HgCdTe  alloys  are  reported.  The  calcu¬ 
lated  values  of  TPA  and  FCA  in  InSb,  InAs,  and  HgCdTe 
alloys  (with  a  =  0.24  and  .\  =  0.34)  are  also  fitted  to  the  simple 
analytical  form  so  they  can  be  used  easily  in  light  propaga¬ 
tion  modeling.3  6'-4  These  values  are  used  in  obtaining  a  so¬ 
lution  to  the  rate  equation  that  describes  the  light  propagation 
through  InAs.  The  results  are  then  compared  with  the  infor¬ 
mation  found  in  the  experiments. 
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THEORY 

Second-order  perturbation  theory  is  used  to  calculate  the 
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TPA  and  FCA  coefficients  -  and  the  Auger  and  radiative  re- 
combination  lifetimes.-  The  dipole  matrix  elements  re¬ 
quired  to  study  light  interaction  with  matter  and  the  scatter¬ 
ing  matrix  elements  required  in  lifetime  calculations  are 
calculated  from  the  underlying  band  structures.  Care  is  given 
to  describe  the  occupancy  of  initial  and  final  states  with  ap¬ 
propriate  Fermi-Dirac  distributions  and  to  include  an  accu¬ 
rate  sum  over  crystal  wave  vectors  spanning  the  entire  Bril- 
louin  zone  (BZ). 

Band  structures 

The  calculations  reported  here  are  based  on  a  hybrid 
pseudopotential  tight-binding  Hamiltonian.  In  this  formal¬ 
ism,  the  tight-binding  Hamiltonian  contains  two  parts.  The 
dominant  part  is  derived  from  the  empirical  pseudopotential 
and  includes  all  long-range  interactions.  A  short-range  em¬ 
pirical  Hamiltonian  is  added  to  fine  tune  the  calculated  band 
structure  so  that  it  agrees  with  the  measured  energy  gaps  and 
effective  masses  at  various  symmetry  points  in  the  BZ.  In  the 
sp 3  basis,  with  the  spin-orbit  coupling  included,  the  Hamil¬ 
tonian  takes  the  16  X  16  matrix  form.  In  the  case  of  alloys, 
we  find  that  the  correction  to  bands  from  the  alloy  disorder  is 
small  near  the  band  edge,  and  hence,  the  virtual  crystal  ap¬ 
proximation  can  be  used  to  obtain  the  alloy  band  structures. 
The  resulting  Hamiltonian  matrix  is  diagonalized  to  yield  the 
eigenvalues  (energies)  and  eigenvectors  (wave  functions). 

Since  high  intensity  light  absorbed  in  a  material  can 
raise  the  material  temperature  substantially  above  room  tem¬ 
perature,  band  structures  over  a  wide  range  of  temperatures 
from  77  K  up  to  —500  K  need  to  be  used.  Although  the 
temperature  (7)  dependence  of  the  matrix  elements  can  be 
included-  for  the  calculation  of  the  band  structure,  those 
calculations  are  extremely  time-consuming,  particularly 
when  the  integration  over  the  full  BZ  is  required  for  lifetime 
calculations.  Here,  as  an  alternative,  the  effect  of  T  on  the 
Hg!_tCdtTe  alloy  band  structure  is  included  indirectly  as  fol¬ 
lows.  First,  for  a  given  T  and  Cd  concentration  x,  the  band 
gap  Eg  in  eV  is  determined  from  the  widely  used 
relationship30  -0.302+ 1.93x+5.35  X  10-4  (1 -2x)*7-0.81x2 
+  0.832x3.  We  then  adjust  the  value  of  x  in  our  zero- 
temperature  Hamiltonian  to  x'  to  get  the  band  gap  of  Eg.  The 
zero-temperature  Hamiltonian  with  the  value  of  x'  is  then 
used  in  further  calculations.  In  the  case  of  InAs  and  InSb 
compounds,  the  parameters  of  the  short-range  empirical 
Hamiltonian  are  adjusted  at  each  temperature  to  yield  the 
band  gap  in  eV  given  by  0.415-2.76 X  10_472/(7’+83)  and 
0.24-6  X  10-472/(7+500),  respectively.31,32  The  calculated 
band  gap  energy  as  a  function  of  T  is  shown  as  a  dashed  line 
in  Fig.  1  for  (a)  InAs,  (b)  Hg066Cd0  34Te,  (c)  InSb,  and  (d) 
Hgo.76Cdo.24Te.  We  see  that  the  band  gap  decreases  with  T  in 
InAs  and  InSb,  but  increases  in  HgCdTe  alloys. 

Fermi  level 

The  study  of  electron-photon  interaction  requires  the 
knowledge  of  the  electron  density  and  its  distribution.  The 
Boltzmann  distribution  is  often  used.  However,  at  high  tem¬ 


peratures  and  for  small  gap  materials,  the  bands  are  nearly  or 
fully  degenerate  and  the  Fermi-Dirac  (FD)  distribution  func¬ 
tion, 

/(£,£>)  =  [1  +  e(E-EfVkBTj-\  (i) 

needs  to  be  used,  where  kB  is  the  Boltzmann  constant.  The 
equilibrium  Fermi  level  (FL)  or  Fermi  energy  £)  in  Eq.  (1)  is 
obtained  from  the  charge  neutrality  condition  that  the  total 
positive  charges  (the  sum  of  equilibrium  hole  density  p0  and 
ionized  donor  density  /Vj)  equals  the  total  negative  charges 
(the  sum  of  equilibrium  electron  density  nQ  and  ionized  ac¬ 
ceptor  density  /V)j).”  The  quantities  n0,  p0 ,  /Vj,  and  /Vj,  in 
turn,  depend  on  the  Fermi  energy  as  follows: 

=  Nd[  1  +  ge(h:rEr>Vkn'y , 

(2) 

N~a  =  Na[  1  +  1 , 

/•oo 

«0  =  2  pc(E-Ec)f{E,Ef)dE, 

Jec 

fE„ 

P  0  =  2  Pv(E-Ev)f(E,Ef)dE,  (3) 

J-  00 

S(Ekc  v-E), 

where  Q  is  the  crystal  volume,  the  degeneracy  factor  g  is  2 
for  CB  and  4  for  VB,  and  Ec  and  Ev  are  the  band  energies  at 
CB  minimum  and  VB  maximum,  respectively.  In  the  special 
case  where  the  material  is  intrinsic  ( ND=NA  =  0 ),  the  FL  is 
denoted  by  fij  and  the  corresponding  carrier  density  (n0 
=Po)  is  denoted  by  nr  The  calculated  E'j-  is  shown  by  the 
dotted  line  in  Fig.  1  for  (a)  InAs,  (b)  Hg0  .66Cd0.34Te,  (c) 
InSb,  and  (d)  Hg0  76Cd0  24Te.  All  energies  are  measured  with 
respect  to  the  top  of  the  VB  at  that  temperature.  At  low 
temperature,  the  E'j  is  near  the  midgap.  Owing  to  the  differ¬ 
ence  in  the  effective  masses  of  electrons  and  holes,  the  FL 
moves  closer  to  the  CB  edge  (shown  as  a  dashed  line)  with 
an  increase  in  T.  In  the  case  of  InSb,  we  see  that  the  FL  lies 
above  the  bottom  of  the  CB  at  high  temperature,  making  the 
material  fully  degenerate.  In  HgCdTe  alloys,  the  band  gap 
increases  with  T  and  hence  the  intrinsic  FL  continues  to  lie 
below  the  CB  for  all  the  temperatures  considered  here.  The 
corresponding  intrinsic  carrier  density  n,  obtained  from  Eq. 
(3)  is  plotted  in  Fig.  2  for  all  four  materials.  The  intrinsic 
carrier  density  increases  with  T  but  decreases  with  Eg  and 
varies  approximately  as  exp(-Eg/2kBT).  This  behavior  is 
clearly  seen  in  Fig.  2.  Recently,  the  variable  Hall  studies,  in 
combination  with  multiple  carrier  analysis,  have  demon¬ 
strated  the  ability  to  yield  an  accurate  measurement  of  the 
carrier  density  in  doped  samples.34  36  The  carrier  densities 
calculated  using  our  band  structures  were  found  to  be  in 
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excellent  agreement  with  these  measurements. “ 

Under  high  intensity  light  illumination,  photoexcited 
electron  density  (Ah)  and  hole  density  (A p)  are  created,  with 
An  =  A p,  which  can  be  far  larger  than  the  equilibrium  carrier 
density  n0  and  p0.  Consequently,  the  TPA,  FCA  coefficients, 
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FIG.  1.  Variation  of  energy  gap  (dashed  line)  and  intrinsic  Fermi  energy  (dotted  line)  with  temperature.  Also  shown  are  electron  quasi-Fermi  energy  (thin  solid 
line)  and  hole  quasi-Fermi  energy  (thick  solid  line)  for  various  photoexcited  carrier  densities  (A«=Ap). 


and  lifetime  calculations  have  to  be  carried  out  not  only  un¬ 
der  degenerate  conditions  but  also  under  nonequilibrium 
condition,  in  which  the  electrons  in  the  CB  and  the  holes  in 
the  VB  have  different  Fermi  energies.  This  quasi-Fermi  en¬ 
ergy  for  electrons  (Ej)  and  holes  (Ej)  is  determined  from  the 
following  conditions: 

r  oo 

n0  +  An  =  2  pc(E  -  Ec)f(E,EfidE, 

(4) 

rEv 

Po  +  Ap  =  2  pc(E  -  Ev)f(E,E})dE. 

J  _oo 

Assuming  that  the  material  is  intrinsic,  we  have  calcu¬ 
lated  the  quasi-FLs  for  various  values  of  the  photoexcited 


carrier  density  using  Eq.  (4).  Also  shown  in  Fig.  1  are  cal¬ 
culated  electron  (thin  solid  line)  and  hole  (thick  solid  line) 
Fermi  energies  as  functions  of  the  lattice  temperature  and 
excited  carrier  density.  When  the  photoexcited  carrier  density 
(An=Ap)  becomes  appreciably  larger  than  the  intrinsic  car¬ 
rier  density  «,■  at  that  T,  the  quasi-FL  for  electrons  (holes)  is 
closer  to  the  CB  (VB).  However,  when  the  n,  is  larger  than 
An  or  A p,  the  quasi-FLs  merge  with  the  intrinsic  FL  as 
shown  in  Fig.  1.  The  charge  carrier  distribution  described  by 
the  FD  function  differs  considerably  from  that  described  by 
the  Maxwell-Boltzmann  function  when  the  quasi-FLs  are 
closer  to  the  band  edge.  Consequently,  we  see  from  Fig.  1 
that  the  quasi-FLs  play  a  major  role  even  in  systems  with 
moderate  excited  carrier  densities  at  low  T  and  in  systems 
with  large  densities  at  high  T,  in  general. 
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FIG.  2.  Variation  of  intrinsic  carrier  density  with  temperature. 


Two-photon  absorption 


Electrons  in  the  heavy-hole  (HH)  and  light-hole  (LH) 
valence  bands  absorb  two  photons  to  reach  the  CB  states.  In 
the  second-order  perturbation  theory,  the  absorption  rate  is 
proportional  to  the  occupation  of  the  initial  state  in  the  VB, 
availability  of  the  final  state  in  the  CB,  and  the  dipole  matrix 
element  between  the  VB  and  CB  states  at  a  given  crystal 
momentum.  Explicitly,  the  expression  for  TPA  coefficient  [3 
is 


P-- 


1  /  2m\ 


7 TC 


i2  [1 


m  /  (w«)3"kl 

X \MkJ2f(Eka,Ef)  S(Eka  -Ek  +  2hoj), 


(5) 


M 


aft ' 


=  22 

P'  jj' 


<^le/  •  pl^'X^'kj  •  pl^g) 

Ek-Ek+hco 


(6) 


where  E^icffy  is  the  electron  energy  (wave  function)  in  band 
/ 3  at  the  wave  vector  k  in  the  BZ,  e  and  m  are  the  free 
electron  charge  and  mass,  n  is  the  real  part  of  the  refractive 
index,  p  is  the  electric  dipole  vector,  ej  is  the  unit  vector  in 
the  polarization  direction  of  light,  and  ha>  is  the  photon  en¬ 
ergy.  The  distribution  function  for  the  electrons  /( Ek_,)  given 
by  Eq.  (1)  uses  appropriate  quasi-Fermi  energies  EJ  and  Ej 
depending  on  whether  f3  denotes  a  valence  or  a  conduction 
band.  In  Eq.  (5),  a  and  [3  denote  respectively  the  initial  and 
final  states  in  any  occupied  VB  and  any  unoccupied  CB  that 
satisfies  the  energy  conservation  condition  specified  by  the  S 
function.  The  sum  over  k  is  performed  over  the  entire  BZ.  A 
careful  analysis  shows  that  the  intermediate  state  f3'  can  be 
any  of  the  unoccupied  CB  state  and  one  VB  state  from  which 
the  electron  excitation  begins.  In  fact,  the  leading  term  of  Eq. 
(6)  arises  when  /?'  is  the  same  as  a  In  that  case,  the  denomi¬ 
nator  of  Eq.  (6)  is  simply  ha>,  ( cffy e  j  ■  •  p </k)  is  nearly  a  con¬ 
stant  for  k  close  to  the  BZ  center,  and  {(jk)\C.y  •  p|  (fy  is  pro¬ 
portional  to  k.  After  substituting  for  Majg  in  Eq.  (5)  under  the 
assumption  that  VB  is  occupied,  CB  is  unoccupied,  and  the 
bands  are  described  in  parabolic  approximation,  it  is  straight¬ 
forward  to  show  that  [3  is  proportional  to  (2hoj- E,,)212  / a>5. 


This  functional  dependence  of  f3  is  being  widely  used  in  the 

literature.6’9’11’12'14-17’20 

We  used  full  band  structures  to  accurately  calculate  the 
FD  distribution  function,  band  energy,  wave  functions,  and 
optical  matrix  elements  to  obtain  /3  as  a  function  of  the  pho¬ 
ton  wavelength  X  and  T.  The  calculated  f3  at  7  =  300  and 
77  K  are  shown  in  Fig.  3  for  (a)  InAs,  (b)  Hg0  .66Cdo.34Te,  (c) 
InSb,  and  (d)  Hg076Cd024Te.  In  all  cases  considered  here,  f3 
increases  with  X,  reaches  a  maximum,  and  decreases  to  zero. 
The  leading  term  in  Eq.  (6)  varies  as  <w_1  and,  consequently, 
/3  in  Eq.  (5)  varies  as  a>~5  (or  X5)  (3  in  Fig.  3  increases 
initially  with  the  wavelength  X.  However  f3  is  also  propor¬ 
tional  to  the  density  of  states  at  the  excitation  energy 
(2hc/k- E,,)  in  the  conduction  band.  The  density  of  states 
decreases  when  X  is  increased.  The  competition  between 
these  two  terms  results  in  a  maximum  in  f3.  When  X  is  longer 
than  twice  the  cutoff  wavelength  \g  (with  \  =hd Eg),  elec¬ 
trons  in  the  VB  cannot  reach  the  CB  and  (3  drops  to  zero.  The 
approximate  value  of  X  where  the  f3  reaches  a  maximum  can 
be  estimated  from  Eqs.  (5)  and  (6).  In  the  approximation 
where  the  dipole  matrix  element  is  proportional  to  k,  the 
initial  state  is  always  occupied,  the  final  state  is  always 
empty,  and  the  bands  can  be  described  by  an  effective  mass, 
1 3  from  Eqs.  (5)  and  (6)  is  proportional  to  k5(2hc/\—Eg)3/2. 
The  maximum  of  this  function  occurs  when  X  is  1 .4X,,.  In 
general,  the  band  structure  is  not  parabolic.  In  our  full  cal¬ 
culations,  we  find  that  TPA  reaches  its  largest  value  when  X 
is  approximately  1.5Xg.  This  empirical  relationship  holds 
well  for  all  four  materials  and  two  temperatures  considered 
here. 

As  the  one-photon  absorption  will  start  to  dominate  at 
wavelengths  shorter  than  the  cutoff  wavelength  for  the  given 
T,  the  values  of  (3  are  not  calculated  in  that  spectral  region. 
In  InAs  and  InSb,  the  band  gap  decreases  when  T  is  in¬ 
creased  and  hence  the  density  of  states  (DOS)  at  the  excita¬ 
tion  energy  (2 hc/k-Eg)  increases,  resulting  in  a  larger  (3  at 
higher  T.  For  similar  reasons,  f3  in  HgCdTe  decreases  when 
T  is  increased.  Although  the  band  gaps  of  InAs  and 
Hgo.eeCdo.  34Te  at  room  temperature  are  equal,  notice  that  [3 
in  the  HgCdTe  alloy  is  considerably  higher  than  in  InAs 
owing  to  the  difference  in  their  density  of  states. 

For  an  intrinsic  material  at  low  temperature,  the  initial 
state  is  full  and  the  final  state  is  almost  empty.  However,  as 
discussed  in  the  previous  section,  the  occupancy  or  availabil¬ 
ity  of  the  states  depends  not  only  on  T  but  also  on  the  excess 
carrier  density.  Consequently,  the  calculated  (3  has  a  depen¬ 
dence  on  the  photoexcited  carrier  density,  as  shown  in  Fig.  4 
at  (a)  300  K  and  (b)  77  K.  At  300  K,  the  intrinsic  carrier 
density  in  InAs  and  Hg066Cd034Te  is  small  because  of  the 
large  band  gap,  and  the  calculated  f3  shows  only  a  slight 
variation  with  the  photoexcited  carrier  density.  However,  in 
the  small  gap  InSb  and  Hg07()Cd0  24Te,  the  intrinsic  carrier 
density  is  large.  Additional  photoexcited  carriers  deplete  (fill) 
the  near  band  edge  VB  (CB)  states.  Hence,  the  occupancy 
(availability)  of  the  initial  (final)  state  is  different  from  1 
(zero),  resulting  in  a  decreased  j3.  We  see  that  the  calculated 
/ 3  starts  to  decrease  for  the  photoexcited  carrier  densities  in 
excess  of  mid  1016  cm-3.  This  effect,  commonly  known  as 
the  Moss-Burstein  effect,  is  dominant  for  small  band  gap 
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FIG.  3.  Two-photon  absorption  coefficient  as  a  function  of  wavelength  and  temperature. 


materials  and  at  low  T.  The  fi  calculated  at  77  K  shows  an 
even  stronger  carrier  density  dependence,  as  shown  in  Fig. 
4(b). 

It  is  instructive  to  compare  the  calculated  /3  to  that  ob¬ 
tained  by  most  commonly  used,  but  approximate,  methods. 
Qualitatively,  the  wavelength  dependence  shown  in  Fig.  4 
agrees  well  with  that  predicted  by  Wherrett’s  or  Van  Stry- 
land’s  expression.8'11  Our  calculated  values  and  those  ob¬ 
tained  with  Van  Stryland’s  expression  are  given  in  Table  I  for 
two  temperatures  and  wavelengths.  We  see  that  our  values  of 
fi  calculated  at  4.8  pm  and  at  low  carrier  density  are  smaller 
than  those  obtained  with  Van  Stryland’s  expression  by  a  fac¬ 
tor  of  2.5  for  Hg0  66Cd0.34Te  and  by  1.5  for  InAs.  Similarly, 
the  values  of  fi  calculated  at  9.6  /xm  and  low  carrier  density 
are  smaller  by  a  factor  of  2.2  for  InSb  and  by  1.2  for 

Hgo.74Cdo.26Te. 

Unlike  InAs  and  InSb,  the  band  gap  of  the  Hg|_lCdvTe 
alloy  can  be  continuously  varied  from  0  to  1.6  eV  by  chang¬ 
ing  Cd  concentration.  Since  this  alloy  system  can  be  used  to 
study  the  propagation  of  light  at  a  number  of  wavelengths, 
the  wavelength-dependent  calculations  of  fi  at  low  carrier 


densities  are  carried  out  as  a  function  of  x  and  are  shown  in 
Fig.  5  for  (a)  T=  300  K  and  (b)  T=ll  K.  The  wavelength 
dependence  is  similar  in  all  alloys  at  both  temperatures. 
Flowever,  as  the  value  of  x  increases,  the  band  gap  increases 
and,  consequently,  fi  decreases.  Similarly,  the  band  gap  in¬ 
crease  with  temperature  also  causes  the  decrease  in  f3.  This  is 
clearly  seen  in  Fig.  5.  Interestingly,  [3  reaches  the  maximum 
when  1.5X.J,  for  these  alloys  and  temperatures,  except  for 
x=0.2  alloy  at  77  K.  The  maximum  occurs  when  \  ~  1 ,2k,,. 
Since  the  band  gap  is  very  small  (83  meV)  for  this  alloy  at 
77  K,  the  conduction  band  near  the  BZ  center  is  mostly  lin¬ 
ear  in  k,  and  not  parabolic.  For  these  bands,  fi  is  proportional 
to  \5(2hcl\-Eg)2  and  the  maximum  of  this  function  is  at 
1.2A.„,  in  agreement  with  our  calculated  value. 


Free  carrier  absorption 

The  photocarriers  excited  by  the  TPA  can  cause  further 
absorption  of  light.  When  a  hole  is  photogenerated  in  the 
heavy-hole  VB,  the  electrons  in  the  light-hole  VB  can  absorb 
one  photon  and  get  excited  to  fill  that  hole.  This  direct  pro- 
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FIG.  4.  Two-photon  absorption  coefficient  as  a  function  of  excited  carrier 
density  An  at  (a)  300  K  and  (b)  77  K.  The  right  axis  in  (b)  corresponds  to 
x=0.24  alloy. 


FIG.  5.  Two-photon  absorption  coefficient  as  a  function  of  wavelength  in 
HgCdTe  alloys  at  (a)  300  K  and  (b)  77  K. 


cess  is  normally  known  as  FCA  by  holes.  The  photogener¬ 
ated  electron  in  the  CB  can  also  absorb  a  photon  through  an 
indirect  process.  To  satisfy  the  momentum  and  energy  con¬ 
servation  conditions,  the  photon  absorption  by  conduction 
electrons  should  be  accompanied  by  a  phonon  absorption  or 
emission.  The  free  carrier  absorption  rate  calculated  from 


TABLE  I.  The  TPA  coefficients  calculated  with  full  bands  are  compared 
with  those  obtained  using  the  analytical  expression. 


Parameters 

X 

(/4m) 

T 

(K) 

/ 3  (analytical) 
(cm/MW) 

/3  (our  work) 
(cm/MW) 

InAs 

4.8 

77 

0.74 

0.31 

Hgo.66Cdo.34Te 

4.8 

77 

2.09 

1.37 

InSb 

9.6 

77 

1.66 

0.81 

Hgo.76Cdo.44Te 

9.6 

77 

15.11 

12.23 

InAs 

4.8 

300 

1.25 

0.49 

Hgo.66Cdo.34Te 

4.8 

300 

1.63 

1.04 

InSb 

9.5 

300 

8.37 

3.30 

HgO.76Cdo.44Te 

9.5 

300 

5.71 

4.68 

second-order  perturbation  theory  is  used  to  obtain  the  FCA 

coefficient  af  and  the  FCA  cross  section  cr.  Although  the 
J  22 
basic  formalism  has  been  described  previously,  "  here  the 

explicit  expressions  for  two  contributions — the  direct  term 

and  the  indirect  term — are  provided.  This  helps  us  separate 

the  contributions  from  electrons  and  holes  for  appropriate 

use  in  the  light  propagation  equations.  In  the  direct  term,  an 

electron  from  a  lower  band  absorbs  a  photon  and  transfers  to 

a  higher  band.  ’  The  initial  and  final  state  wave  vectors 

remain  unchanged,  but  the  final  state  energy  is  larger  than  the 

initial  state  energy  by  hco0.  We  have 


aD-- 


2i re\2/  1  \ 


X 


m  /  \fio)0c/a/lk 


2  f(Eka,E")[] -f(E%Ef)] 


S(Ek-Eka-hu  0). 


(7) 


Since  the  energy  separation  between  the  LH  and  HH  bands 
varies  from  0  to  several  eV,  this  term  is  nonzero  for  hole- 
initiated  FCA.  However,  the  CBs  are  separated  by  energies 

6 


Downloaded  07  Nov  2007  to  134.131.125.49.  Redistribution  subject  to  AIP  license  or  copyright,  see  http://jap.aip.org/jap/copyright.jsp 


113104-7  Krishnamurthy  et  a/. 


J.  Appl.  Phys.  101,  113104  (2007) 


much  larger  than  the  usual  photon  energies  of  interest,  and 
the  direct  excitation  of  electrons  from  the  lowest  CB  to  the 
next  highest  CB  is  not  possible.  Hence,  in  Eq.  (7)  both  a  and 
/3  denote  valence  bands  and  the  FD  function  is  evaluated 
with  EVj.  An  additional  mechanism  to  supply  both  momentum 
and  energy  is  required  for  CB  electrons  to  participate  in  the 
FCA.  The  most  dominant  contribution  usually  arises  from 
longitudinal-phonon  (FO)  assisted  FCA.  Using  the  Frohlich 
Hamiltonian  to  describe  electron-phonon  interaction,  the  ab¬ 
sorption  coefficient  for  this  indirect  process  can  be  written  as 


(27r)3e4fi(oLO 

m2hc(h(o)3 


2  f{El,Ej) 

a  (3k 


X2  [1  -f{Ef\Efi\R^dLEl±huh0  +  hu-Ef(ll 


*%>" 


1  _  1 

Nlo+2  +  2 


k±q 


o'*T 


(8) 


$/2  «<^lej-p|<#;> 


-<<q|erp|^±q» 


2 


where  the  ha>L0  and  /VLO  are  the  FO  phonon  energy  and 
phonon  population,  k0  and  k.,  are  the  dielectric  constant  at 
zero  and  infinite  frequency,  and  the  ±  sign  denotes  phonon 
absorption  or  emission.  For  the  conduction  electrons,  the  ini¬ 
tial,  intermediate,  and  final  states  will  all  be  in  the  lowest 
CB,  whereas  for  the  holes,  the  initial  and  final  states  can  be 
in  either  the  FH  or  HH  bands.  The  FD  functions  are  evalu¬ 
ated  using  appropriate  quasi-Fermi  level  depending  on  the 
bands  that  contribute  to  the  absorption.  The  total  FCA  coef¬ 
ficient  ctj?  is  the  sum  of  the  two  contributions  given  by  Eqs. 
(7)  and  (8),  and  the  FCA  cross  section  for  electrons  (holes) 
cre(;i)  is  the  ratio  of  the  coefficient  to  the  total  electron  (hole) 
density  Ne(hy  Simply, 


af=  aD  +  afD  +  aID, 

(9) 


^e{h) 


af 

^e(h) 


The  calculated  rate  for  FCA  by  electrons,  which  has  only 
the  indirect  term,  is  much  smaller  than  that  for  the  holes.  The 
FCA  coefficients  for  both  electrons  and  holes  were  found  to 
be  proportional  to  the  carrier  density.  Consequently,  cr  is  in¬ 
dependent  of  the  carrier  density.  The  calculated  values  of  cre, 
crh,  and  the  total  cross  section  cr  are  plotted  as  a  function  of 
wavelength  in  Fig.  6  for  intrinsic  InAs  and  Hg0  ,66Cd0.34Te 
alloys.  As  expected  the  electron  contribution  is  much  smaller 
than  the  hole  contribution.  The  total  cross  section  is  small  at 
shorter  wavelengths  because  of  the  unavailability  of  holes  at 
those  energies  in  the  HH  band.  Since  the  hole  density  is 
greater  near  the  zone  center,  the  FCA  cross  section  increases 
with  the  photon  wavelength.  However,  because  of  the  re¬ 
duced  intermixing  of  symmetries  in  both  FH  and  HH  band 
states  near  the  zone  center,  the  dipole  matrix  element  de¬ 
creases  (and  reaches  zero  at  the  center).  Consequently,  the 
cross  section  reduces  with  further  increase  in  wavelength. 


2  4  6  8  10  12  14 

Wavelength  [jjm] 


Wavelength  him] 

FIG.  6.  Wavelength-dependent  free  carrier  absorption  cross  section  at 
300  K. 

The  wavelength  dependence  of  the  FCA  cross  section  in 
InSb  and  Hg066Cd034Te  alloys  (not  shown  in  here)  is  also 
similar. 

When  light  energy  is  absorbed  in  semiconductors,  equal 
densities  of  excited  electrons  and  holes  are  created,  which  in 
turn  release  that  energy  to  the  lattice  in  the  form  of  heat.  To 
model  light  propagation  through  semiconductors,  it  is  there¬ 
fore  necessary  to  know  the  FCA  cross  section  as  a  function 
of  both  carrier  concentration  and  temperature.  The  cr  calcu¬ 
lated  for  four  materials  as  a  function  of  wavelength  at  two 
temperatures  are  shown  in  Fig.  7.  In  the  literature  cr  is  usu¬ 
ally  assumed  to  be  independent  of  temperature  T.  However, 
Fig.  7  shows  that  cr  varies  strongly  with  T — at  long  wave¬ 
lengths,  cr  decreases  with  T,  while  at  short  wavelengths,  it 
increases  with  T.  This  dependence  can  be  explained  by  ex¬ 
amining  the  transition  rate  obtained  using  the  second-order 
perturbation  theory.  The  absorption  coefficient  is  approxi¬ 
mately  the  integrated  product  of  the  hole  density  and  optical 
matrix  element.  Although  this  coefficient  is  divided  by  the 
total  hole  density  to  get  cr,  only  a  fraction  of  the  holes  par- 
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FIG.  7.  Wavelength-dependent  free  carrier  absorption  cross  section  at  77  and  300  K. 


ticipate  in  the  absorption  process  because  of  energy  and  mo¬ 
mentum  conservation  conditions.  Even  when  the  matrix  ele¬ 
ments  do  not  change  with  T,  the  ratio  of  participating  holes 
to  the  total  number  of  holes  has  a  strong  dependence  on  T. 
For  larger  photon  energies  (or  shorter  wavelengths),  the  VB 
states  away  from  the  zone  center  participate  in  absorption 
and  the  hole  density  in  those  states  increases  with  T.  For 
longer  wavelength  photon  absorption,  the  holes  near  the  zone 
center  participate  and  the  hole  density  decreases  with  T  as 
the  holes  are  thermally  excited  to  higher  energies.  Conse¬ 
quently,  cr  decreases  with  T. 

We  have  calculated  the  T  dependence  of  cr  in  InAs  and 
Hgo.66Cd0.  34Te  at  a  wavelength  of  4.8  /xm  and  in  InSb  and 
Hgo.veCdo.  24'fe  at  a  wavelength  of  9.6  /im.  Figure  8  shows 
that  the  cr  at  9.6  jim  decreases  by  a  factor  of  14-16  as  the  T 
is  increased  from  77  to  500  K.  cr  at  4.8  /cm  in 
Hgo.66Cd0.  34Te  is  nearly  a  constant,  but  increases  rapidly 
with  T  in  In  As.  It  is  important  to  note  that  the  T  dependence 
of  cr  is,  in  general,  very  strong  except  at  one  or  two  wave¬ 
lengths  where  the  hole  density  does  not  change  with  T. 

The  wavelength  dependence  of  cr  for  various  x  values  in 
Hg|_vC’dvTe  alloys  are  shown  in  Fig.  9  for  (a)  7  =  300  K  and 


(b)  T=  77  K.  As  x  increases,  the  band  gap  increases  and  the 
hole  density  decreases.  Hence  cr  is  expected  to  decrease 


100  200  300  400  500  600 

T  [K] 


FIG.  8.  Temperature-dependent  free  carrier  absorption  at  4.8  /zm  in  InAs 
and  Hg066Cd034Te  and  at  9.6  /zm  in  InSb  and  Hg076Cdo24Te  at  300  K. 
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Wavelength  [pm] 


Wavelength  [pm] 

FIG.  9.  Wavelength-dependent  free  carrier  absorption  cross  section  as  a 
function  of  Cd  concentration  in  HgCdTe  alloys  at  (a)  300  K  and  (b)  77  K. 


The  lifetimes  limited  by  the  Auger  recombination  (AR)  and 
radiative  recombination  (RR)  are  obtained  as  functions  of 
temperature  and  carrier  densities. 


Parameter  fit 

Although  the  light  propagation  modeling  studies3-6'22"4 
use  only  a  few  parameters — /?,  o\  t,  N(=n0+An),  T,  and 
sample  thickness  d — complications  arise  because  of  their  in¬ 
terdependence.  For  example,  /?  is  a  function  of  both  N  and  T, 
but  N  depends  on  f3  and  T.  Similarly,  r  depends  on  T  and  N, 
which  in  turn  depends  on  T.  Eventually,  /?,  cr,  and  t  will 
determine  the  increase  in  T,  which  will  in  turn  decide  the 
values  of  those  input  parameters.  Hence,  for  a  self-consistent 
solution  to  the  rate  equations,  it  is  preferable  to  have  a  closed 
functional  form  for  the  parameters  /3,  cr,  and  t,  which  is  valid 
in  the  temperature  and  carrier  density  ranges  of  interest.  The 
values  of  (3  and  cr  have  been  calculated  at  a  wavelength  of 
4.8  /im  in  InAs  and  Hg0  .66Cd().  34Te  and  at  9.6  pm  in  InSb 
and  Hg0  V6Cd0  24Te,  for  77  K  <7’<500  K  and  10 14  cm-3 
<  A n<  1018  cm-3.  They  are  then  fitted  to  the  following  ex- 
pressions: 

/3(cm/MW)  =  — - £ - — , 

1  +  expL(ft  -  c)/b\ 

4  4 

b  =  21  bit,  c  =  2  C,T , 

/= 0  i= 0 

n  =  (N  X  1CT18), 


cr(10  16  cm2)  =  2  o-jV, 
i= o 


!=0 


(10) 


(ID 


when  x  is  increased.  However,  as  x  is  increased,  the  spin- 
orbit  coupling  increases  and  the  LH  band  bends  more  away 
from  the  HH  band.  Consequently,  the  energy  conservation 
condition  is  satisfied  in  the  BZ  closer  to  the  zone  center.  This 
increases  the  density  of  participating  holes,  and  hence,  cr 
expected  to  increase  with  x  is  increased.  The  relative 
strengths  of  these  two  competing  mechanisms  depend  on  the 
wavelength  and  temperature.  The  crossover  feature  is  seen 
clearly  more  at  77  K  than  at  300  K. 


Lifetimes 

The  lifetimes  t  of  the  excited  electrons  play  a  vital  role 
in  light  propagation  through  the  material  because  the  excited 
carriers  can  participate  in  further  absorption  of  light  only 
before  they  recombine.  The  full  band  structures  and  accurate 
scattering  matrix  elements  are  used  to  obtain  lifetimes  as 
described  previously."  Those  calculations  are  here  gener¬ 
alized  to  use  the  quasi-Fermi  levels  for  electrons  and  holes. 


rAR(ns)  =  10(*+18,V,  g  =  ^gir,  h  =  ^  h,r,  (12) 

/= 0  /=0 

..  4^  (  1000 y 

A,(1015  cm  3)  =  e",  n  =  ^lni\——  .  (13) 

N  is  the  total  carrier  density,  which  is  the  sum  of  n0  and  A« 
(or  pq  and  A p  for  holes)  in  the  unit  of  cm-3.  The  parameters 
n,  /;,  c,  g,  h,  and  n  are  dimensionless.  The  parameters  fitted 
to  the  calculated  values  are  given  in  Table  II  for  the  four 
materials  considered  here.  The  table  lists  parameters  only  for 
AR  lifetimes.  We  had  previously  calculated  RR  lifetimes  and 

25 

found  them,  in  general,  to  be  very  long  for  these  materials. 
In  addition,  because  of  the  photon  recycling  present  in  the 
thick  samples  studied  here,  the  RR  lifetimes  are  extremely 
long  and  have  a  negligible  effect  on  light  propagation.  The 
fitted  parameters  for  RR  lifetimes  are  therefore  not  listed. 
The  AR  lifetimes  in  Hg0  .76Cd0.  24Te  are  not  calculated  and  not 
listed  as  they  are  not  expected  to  be  substantially  different 
from  those  in  Hg0  ,66Cd0.34Te. 
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TABLE  II.  The  parameters  used  to  fit  the  calculated  /?,  cr,  r,  and  A,  at  the  wavelength  of  4.8  /xm  for  InAs  and 
Hgo.66Cdo.34Te  and  at  9.6  fjL m  for  InSb  and  Hg0  76Cd0  24Te,  in  the  temperature  range  from  75  to  500  K  and  in  the 
excited  carrier  density  range  from  of  1014  to  1018  cm-3. 


Parameters 

InAs 

Hgo.66Cdo.34Te 

InSb 

Hg0.76Cd0.44Te 

0.205  16 

2.1620 

2.5913 

48.074 

a 1 

0.001  529  9 

-0.003  947  30 

0.003  969  1 

-0.455  16 

«2 

-1.859  X10“6 

1.3372  X10“5 

0.000  125  89 

0.002  999  3 

a3 

1.8466  X  10“9 

-3.1219X  10-8 

-4.6698  X  10“7 

-8.3954  X  10“6 

a4 

3.4373  X  10-12 

3.3459  X10-11 

6.6038  X  lO"10 

7.718  X  10“9 

bo 

-0.034  375 

0.117  42 

-0.001  972  8 

0.3932 

bi 

0.002  064  5 

0.000  413  77 

0.000  312  68 

-0.003  377  1 

bi 

-2.1 108  X  10-6 

5.3520  X  10"6 

5.2291  X10“7 

2.0397  X10“5 

63 

-6.5561  X  10-10 

-6.6517  X10“9 

7.3333  X  10“10 

-5.051  X10“8 

bA 

5.6352  X  10“12 

5.4334  X10“12 

3.0668  x  10“12 

4.563  X10“n 

c0 

0.458  59 

1.2756 

-0.045  22 

0.210  89 

Cl 

-0.001  353  4 

-0.000  958  61 

0.0 

0.003  298  8 

c2 

1.0666  X  10“5 

1.0237  X  10“5 

0.0 

-3.8836  X  10“5 

c3 

-1.1906  X  10"8 

-1.8417  X10“8 

0.0 

1.1845  X10“7 

c4 

4.3216  X  10“12 

-7.1237  X10“12 

0.0 

-1.1285  X  10_1° 

<ro 

-2.8164 

-0.955  01 

46.1 

65.296 

07 

0.044  098 

0.166  77 

-0.2929 

-0.571  62 

02 

5.9962  X10“6 

-0.000  813  82 

0.000  826  67 

0.002  110  9 

03 

-2.2463  X  10“7 

1.4824  X  10“6 

-1.1092X  10-6 

-3.4952  X  10“6 

04 

2.3863  X  lO”10 

-9.422  X  10"10 

5.6047  X  10“10 

2. 1284  X  10-9 

go 

40.467 

43.429 

47.74 

43.429 

£i 

-0.027  37 

-0.067  12 

-0.152  05 

-0.067  12 

Si 

-8.597  X10“6 

2.480  X10“4 

0.000  752  8 

2.480  X  10“4 

Si 

2.2386  X10“7 

-3.7801  X  10“7 

-1.6296  X10“6 

-3.7801  X  10“7 

Sa 

-3.122  X10“10 

1.8071  X10“10 

1.241  X10“9 

1.8071  X10"10 

K 

-2.1984 

-2.3581 

-2.5943 

-2.3581 

K 

0.000  518  36 

0.003  266  5 

0.007  922  6 

0.003  266  5 

hi 

5.578  X10"6 

-1.1833  X  10-5 

-4.0231  X10“5 

-1.1833  X  10“5 

hi 

-2.2914  X  10"8 

1.7764  X  10"8 

8.9063  X10“8 

1.7764  X10“8 

h4 

2.5089  X  10“n 

-7.8921  X  10“12 

-6.8868  X  10"11 

-7.8921  X10“12 

nQ 

11.119 

9.3877 

10.305 

8.605 

»i 

-3.6906 

-2.9946 

-2.4768 

-2.0477 

»2 

0.155  58 

0.154  55 

0.092  403 

0.164  55 

«3 

-0.010  771 

-0.010  242 

-0.004  431  9 

-0.010  89 

«4 

0.000  296  32 

0.000  264  68 

8.7616  X10“5 

0.000  280  82 

EXPERIMENTS  AND  RESULTS 

For  samples  that  have  a  thickness  that  is  smaller  as  com¬ 
pared  to  the  collimated  region  of  the  laser  beam  incident  on 
them,  it  can  be  assumed  that  the  beam  remains  collimated 
even  at  high  intensities.  The  equation  describing  the  space 
and  time  evolution  of  the  pump  intensity  7P  and  probe  beam 
intensity  /pr  can  be  written  as 


-^  =  -(3ll-(aeNe  +  (rhNh)Ip, 
dz 

(14) 

d(An)  /S/p  An 
dt  2  hu>  tr 

(15) 

C^f-  =  -  +  <rMl pr- 

dz 

(16) 

The  first  term  in  Eq.  (14)  is  from  TPA  and  the  second  term  is 
from  FCA.  Ne  is  the  total  electron  density  (which  is  the  sum 
of  n0  and  An),  Nh  is  the  total  hole  density  (which  is  the  sum 


of  pQ  and  A p),  and  tr  is  the  total  recombination  time  limited 
by  AR,  RR,  and  Shockley-Read-Hall  (SRFI)  mechanisms. 
For  an  intrinsic  material,  the  FCA  term  can  be  simplified  to 
rr/V/p.  The  change  in  the  excited  carrier  density  is  a  resultant 
of  the  gain  from  the  TPA  and  the  loss  to  the  recombination  as 
described  in  Eq.  (15).  The  probe  intensity  is  too  small  to 
cause  TPA,  and  the  rate  equation  for  /pr  [Eq.  (16)]  contains 
only  a  FCA  term.  When  the  high  intensity  is  absorbed  either 
by  TPA  or  by  FCA,  the  absorbed  energy  is  eventually  trans¬ 
ferred  to  the  lattice  as  heat.  The  consequent  change  in  the 
lattice  temperature  can  be  obtained  from  the  rate  equation: 

dT  fill  +  (aeNe  +  (rhNh)I?  ^ 

dt  Cy  Cy 

where  cv  is  the  specific  heat  of  the  material.  The  indirect 
term  in  FCA  should  take  into  account  the  phonon-emission- 
induced  increase  and  phonon-absorption-induced  decrease  in 
T.  The  modified  cross  section  cr  includes  the  difference  in 
the  emission  and  absorption  contributions  in  the  indirect 
term  of  the  FCA,  i.e.,  cr^h)  =  (aD+a'lI)-aiD)/Ne^hy  However, 
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FIG.  10.  Comparison  of  normalized  pump  transmission  coefficient  mea¬ 
sured  (data)  at  300  K  in  1015  //-doped  InAs  with  that  calculated  (solid  lines) 
for  four  thicknesses. 

the  direct  term  in  the  FCA  is  considerably  larger  than  the 
indirect  term;  the  cr  is  nearly  the  same  as  cr.  The  parameters 
n0,  f3,  cr,  and  t  depend  on  T.  In  addition,  (3  and  tr  depend  on 
carrier  density,  which  in  turn  depend  on  T.  Hence  Eqs. 
(14)— (17)  will  have  to  be  simultaneously  solved  to  find  the  7P 
and  7pr  at  a  time  1  and  a  position  z  in  the  sample. 

We  used  the  parameters  for  InAs  listed  in  Table  II  and 
solved  Eqs.  (14)— (17)  to  explain  the  results  obtained  in  our 
pump-probe  experiment.  Since  the  details  are  published 
elsewhere/  we  discuss  the  results  only  briefly  here  to  pro¬ 
vide  continuity.  We  set  up  a  nonlinear  energy  transmission 
and  a  pump-probe  experiment  to  determine  the  validity  of 
the  theory  presented.  In  both  experiments,  a  laser  beam  hav¬ 
ing  a  wavelength  of  4.8  /mi.  obtained  by  frequency  doubling 
a  transverse  excited  atmospheric  (TEA)  CCE  laser  to  provide 
128  ns  (full  width  at  e~l  of  the  maximum)  duration  pulses 
with  energies  up  to  10  mJ,  was  focused  on  InAs  samples. 

In  the  first  experiment,  the  amount  of  energy  incident  on 
a  set  of  samples  of  various  thicknesses  was  varied  and  the 
ratio  of  the  energy  transmitted  through  the  sample  to  the 
incident  energy  was  plotted  as  a  function  of  the  incident 
intensity.  Care  was  taken  to  ensure  that  the  total  energy 
transmitted  by  the  samples  was  collected  in  the  pyroelectric 
energy  detectors  by  using  a  lens  to  focus  the  transmitted 
beam.  The  spatial  distribution  of  the  transmitted  beam  at 
various  intensities  was  measured  at  the  lens  aperture  and  the 
lens  of  appropriate  size  was  chosen  so  that  the  entire  trans¬ 
mitted  energy  was  collected.  First,  nonlinear  absorption  ex¬ 
periments  were  performed  with  only  the  pump  present.  The 
measured  transmitted  intensity  as  a  function  of  incident  in¬ 
tensity  and  sample  thickness  is  shown  in  Fig.  10.  As  the 
intensity  increased,  more  carriers  are  generated  and  the  ab¬ 
sorption  increases,  resulting  in  a  decreased  transmission. 
Similarly,  as  the  thickness  is  increased,  light  is  absorbed 
more.  These  trends  are  clearly  observed.  Also  in  Fig.  10,  we 
have  plotted  the  calculated  transmitted  intensity  obtained 
from  solving  the  rate  equations  [Eqs.  (14)— (17)]  with  the 
values  given  in  Table  II.  In  addition  to  the  Auger  recombi¬ 
nation,  we  have  included  lifetimes  limited  by  the  SRH 


0  100  200  300  400  500x10'® 

Time  [s] 

FIG.  1 1 .  Comparison  of  normalized  probe  transmission  coefficient  mea¬ 
sured  (data)  at  300  K  in  1015  //-doped  InAs  with  that  calculated  (solid  lines) 
for  three  beam  energies. 

mechanism.  The  SRH  scattering  is  not  intrinsic  to  the  mate¬ 
rial  and  depends  on  external  variables  such  as  the  growth 
temperature  and  pressure.  Our  detailed  calculations  carried 
out  previously40  indicate  a  nearly  T-independent  SRH  life¬ 
time  of  200  ns  in  a  2  X  1016  cm-3  //-doped  InAs.  Although 
the  InAs  sample  studied  here  is  only  1015  cm-3  n  doped,  the 
excess  carriers  generated  by  the  pump  is  of  the  order  of 
(1-2)  X  1016  cm-3  and  justifies  the  use  of  the  above  calcu¬ 
lated  SRH  lifetimes.  We  further  assumed  that  the  SRH  life¬ 
time  decreases  linearly  with  excess  carrier  density.  We  see 
that  the  calculated  transmission  intensity  agrees  very  well 
with  measured  values  for  all  thicknesses  studied. 

In  a  second  set  of  experiments,  a  weak  continuous  wave 
probe  beam  (at  a  wavelength  of  5.3  /xm)  was  made  to  be 
incident  on  the  sample  along  with  the  pump  beam  (at 
4.8  /xm )  and  the  temporal  behavior  of  the  probe  beam  was 
monitored  as  a  function  of  the  pump  beam  intensity.  As  the 
pump  intensity  increased,  the  probe  beam  power  dropped 
because  of  the  nonlinear  absorption  in  the  sample.  However, 
at  the  end  of  the  duration  of  the  pump  beam,  i.e.,  after  about 
100  ns,  the  probe  beam  power  recovered  because  of  the  de¬ 
cay  of  the  free  carriers,  and  the  rate  of  recovery  provided  a 
measure  of  the  carrier  lifetime.  In  Fig.  11,  the  experimental 
data  for  the  time  dependence  of  the  probe  beam  are  presented 
for  three  values  of  the  pump  beam  intensity,  along  with  the 
theoretical  calculations.  It  is  seen  that  with  only  the  SRH 
values  used  in  an  approximation,  the  time-dependent  trans¬ 
mitted  probe  intensity  calculated  for  the  three  input  energies 
explains  the  observed  variation  very  well. 

CONCLUSIONS 

We  have  used  a  full  band  structure,  obtained  from  a  hy¬ 
brid  pseudopotential  tight-binding  Hamiltonian,  to  calculate 
the  temperature  and  carrier  density  dependent  two-photon 
absorption  coefficients,  free  carrier  absorption  cross  section, 
and  Auger  lifetimes  in  InAs,  InSb,  and  HgCdTe  alloys.  The 
accuracy  of  the  calculations  is  enhanced  with  the  use  of  the 
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quasi-Fermi  levels,  intrinsic  carrier  densities,  detailed  optical 
and  scattering  matrix  elements,  and  complete  integration 
over  the  entire  BZ.  The  calculated  parameters  are  fitted  to 
closed-form  functions  and  used  in  the  rate  equations  to  ob¬ 
tain  the  transmitted  intensity  as  a  function  of  time,  incident 
intensity,  and  sample  thickness.  The  nonlinear  transmission 
calculated,  with  only  the  SRH  lifetimes  as  an  adjustable  pa¬ 
rameter,  agrees  very  well  the  measured  values.  We  showed 
that  with  accurate  evaluation  of  all  parameters,  the  modeling 
can  provide  detailed  guidance  in  understanding  the  underly¬ 
ing  mechanisms  and  in  choosing  materials  for  light  propaga¬ 
tion  applications.  The  closed-form  expressions  for  the  funda¬ 
mental  parameters  given  here  can  be  used  in  other  accurate 
light  propagation  models. 
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